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Constructing a discrete model like a cellular automaton is a powerful method for understanding various 
dynamical systems. However, the relationship between the discrete model and its continuous analogue is, in 
general, nontrivial. As a quantum-mechanical cellular automaton, a discrete- time quantum walk is defined 
to include various quantum dynamical behavior. Here we generalize a discrete-time quantum walk on a 
line into the feed-forward quantum coin model, which depends on the coin state of the previous step. We 
show that our proposed model has an anomalous slow diffusion characterized by the porous-medium 
equation, while the conventional discrete-time quantum walk model shows ballistic transport. 

Cellular automata - discrete models that follow a set of rules 1 - have been analyzed in various dynamical 
systems in physics, as well as in computational models and theoretical biology; well-known examples 
include crystal growth and the Belousov-Zhabotinsky reaction. To simulate quantum mechanical phe- 
nomena, Feynman 2 proposed a quantum cellular automaton (the Feynman checkerboard). This model, defined 
in the general case by Meyer 3 , is known as the discrete- time quantum walk (DTQW). Since the DTQW on a graph 
is a model of a universal quantum computation 4 ' 5 , it is of great utility, especially in quantum information 6 9 . 
Furthermore, the DTQW has been demonstrated experimentally in various physical systems 10 " 24 to reveal 
quantum nature under dynamical systems. 

As the cellular automaton can be mapped to various differential equations by taking the continuous limit, some 
DTQW models can be mapped to the Dirac equation 25 " 27 , the spatially discretized Schrodinger equation 28,29 , the 
Klein-Gordon equation 27 ' 30 , or various other differential equations 31 ' 32 . These equations have ballistic transport 
properties, which are reflected mathematically in the one- dimensional (ID) DTQW with a time- and spatial- 
independent coin operator, i.e. a ID homogeneous DTQW 33 . We consider here the ID DTQW model. Physically, 
the standard deviation of the homogeneous DTQW is o(t) ~ t, whereas the unbiased classical random walk has a 
standard deviation of a(t) ~ yft. 

In the homogeneous DTQW, the time evolution of a quantum particle (walker) is given by a unitary operator U 
defined on the composite Hilbert space TL S ®7Y C , where Tt s '■ = span{ | j) is the walker Hilbert space, and Tt c 
is the two-dimensional coin Hilbert space. For a unitary operator U, the quantum state evolves in each time step t 
by 



with 



\^ t+1 ) = U\ x ¥ t ) 



j = — GO \ 7 / 



(1) 



(2) 



where the upper aj (lower frj) component corresponds to the left (right) coin state at the j-th site at time step t. As 
an example, the time evolution of the DTQW is given by 



a)t\ = cos 0 Uj — sin 6 bj, 
b'j+l =sin 0 a) + cos Ob). 



(3) 
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The j-th site probability at time step t is given by Pj = aj + 
and y^ 00 p f = i is satisfied for each time step t. 

Z jj = — oo J x 

As a generalization of Eq. (3), we define a DTQW with a 
forward quantum coin described by 



a t+1 - 
a j-i " 



\fy V P 



b t+l - 



with the site- dependent rate function 

«/=k-i|+»'k+i|. 



I |2 



feed- 



(4) 



(5) 



which incorporates the nearest-neighbor interactions. Since this 
quantum coin depends on the probability distribution of the coin 
states on the nearest- neighbor sites at the previous step, this model is 
called & feed- for ward DTQW. It is remarked that the feed-forward 
DTQW is one of the nonlinear DTQW models. Note that if we set the 
rate function tog = cos 0, which is time and site independent, then 
the model in Eq. (4) reduces to the homogeneous model in Eq. (3). 
We will show that our proposed feed-forward DTQW is experiment- 
ally feasible. Furthermore, we will show that this model shows the 
anomalous diffusion as introduced below. 

One of the famous anomalous diffusion equations is the porous 
medium equation (PME) 34 , defined by 



dt 



p(x,t) 



dx 2 



(6) 



where the real parameter m > 1 characterizes the degree of porosity 
of the porous medium. It is known that the PME can be derived from 
three physical equations for the density p, pressure p, and velocity v 
of the gas flow: the equation of continuity, dpldt + V * (pv) = 0; 
Darcy's law, v oc — Vp; and the equation of state for a polytropic gas, 
p oc p v , where v is the polytropic exponent and m = v + 1. One of the 
peculiar features of the PME is the so-called finite propagation, which 
implies the appearance of a free boundary separating the positive 
region (p > 0) from the empty region (p = 0). 

A well-known solution of the PME is the Barenblatt-Pattle (BP) 
one 35 ; it is self-similar, and its total mass is conserved during evolu- 
tion. The evolutionary behavior of the BP solution was recently 
studied in the context of generalized entropies and information geo- 
metry 36 . The BP solution can also be expressed by Tsallis' one-real- 
parameter (q) generalization of a Gaussian function, i.e., the 
^-Gaussian 37 . In the case of ID space, the BP solution is 



Pq(x,t)c 



Eex p fl 



(7) 



with q = 2 — m. Here, a 2 q (t) is a positive parameter that characterizes 
the width of the g- Gaussian at time t and is similar to the variance 
a q= 1 (t) in a standard Gaussian. In other words, the parameter a q (t) 
characterizes the spread of the ^-Gaussian distribution 38 ' 39 ; 

a q (t)cct^, (8) 

which reduces to a q = i (t) oc \ft in the limit of q — > 1. Note that in the 
same limit, the ^-Gaussian reduces to the standard Gaussian, 

exp(-x 2 / ( r^ 1 (0), 

equation dpldt = d 2 pldx 2 . 

In this paper, we analyze a specific feed-forward DTQW with an 
experimental proposal using the polarized state and optical mode. 
We show numerically that the probability distributions of the 
feed-forward DTQW model have anomalous diffusion characterized 



and the PME reduces to the standard heat 



by c q = 05 (t) ~ f A . These dynamics are consistent with the time 
evolution of the self- similar solution 35 of the PME, which is 
known to describe well the anomalous diffusion of an isotropic 
gas through a porous medium. Furthermore, we show analytically 
that the interference terms in our model help the speedup of the 
associated Markovian model but does not help the quadratic 
speedup like the homogeneous DTQW does 40 . Note that although 
anomalous diffusion was found numerically in a nonlinear 
model 41 , an aperiodic time -dependent coin model 42 , and the his- 
tory-dependent coin 43 from the time dependence of the variance 
G q =i(t), the partial differential equation (PDE) corresponding to 
their models have not derived due to the lack of the numerical 
step (about 100 step). Therefore, we have not yet revealed the 
origin of the anomalous diffusion in the DTQW. 

Results 

Experimental proposal of feed-forward DTQW. We propose an 
optical implementation of the feed-forward DTQW. In the simple 
optical implementation of the homogeneous DTQW, the walker 
space uses the spatial mode and the coin space does the polarized 
state. The shift uses the polarized beam splitter and the quantum coin 
uses the quarter-wave, half-wave, and quarter-wave plates, which can 
arbitrarily rotate the polarized state in the Poincare sphere. This was 
experimentally done in Refs. 10-12,16-22. 

Let us construct the feed-forward system of the quantum coin. The 
detectors put at each path to evaluate the probability distribution of 
I 2 I I 2 

the coin state aj and Wa . Since our proposed quantum coin 

depends on |aj and we can calculate the coin operator at the 

jth site. According to the Jones calculation 44 to satisfy Eq. (4), we 
control the angels of the quarter- wave, half-wave, and quarter-wave 
plates for each path. This can be taken as the quantum coin operator 
with the feed-forward. This is depicted in Fig. 1. In what follows, 
we consider the long time time evolution of the feed-forward 
DTQW. 

Numerical results of feed-forward DTQW with anomalous 
diffusion. To study the time evolution of the feed-forward DTQW 
model, the initial state should have nonzero coin states at the nearest- 
neighbor sites. This can be easily understood by considering the 
following example. Let us take (ajj, b° Q ) as the only non-zero initial 
state. In this case, the rate is g® = 0, because there is no neighboring 
state. From the map in Eq. (4), we see that the nonzero states at t = 1 
are a 1 _ l = — b° Q and b\ = ajj. This gives g l _ l =g\= 0, and we see that 
the only nonzero state is (a^bfy = ( — b\,a l _^) = ( — a^, — b^j at t = 
2. This state at t = 2 only differs in sign (or phase) from the initial 
state. Thus if the initial state is concentrated at a single site, no 
spreading occurs; the state only oscillates around the initial site. 

Figure 2 (A) shows a typical probability distribution of the feed- 
forward DTQW after a long-time evolution. See the Supplemen- 
tary Movie for more details. The initial state was set as 
(a° 0 ,b° 0 ) = (a\,b\) = (1/2 ,//2). We note that the probability distri- 
bution diffuses very slowly and does not approach a Gaussian. 
These features are often observed in anomalous diffusion. It is also 
remarked that such behavior has not yet seen in DTQWs with the 
position- dependent coin 45 " 48 , which show the localization property. 

We performed long-time numerical simulations of the feed-for- 
ward DTQW model [Eq. (4)] for up to t ~ 10 8 steps. To study the 
asymptotic behavior, we take running averages of the numerical 
solutions to reduce the influence of multiple spikes. The averaged 
data were fitted with the ^-Gaussian of Eq. (7) to determine the 
corresponding ^-generalized standard deviation a q (t), as shown in 
Fig. 2 (B). We note that the averaged data at each time step are well 
fitted by the Gaussian with q = 0.5. 
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Figure 1 | Optical implementation of the feed-forward DTQW model. Figure shows our experimental proposal of our model. From the intensity of the 
detectors for each path, the polarizers should be changed. This can be taken as the feed-forward quantum coin. 



The long-time evolution of cr q (t), plotted in Fig. 2 (C), reveals that terms. We obtain the following map for both coin distributions 

the time evolution of the feed-forward DTQW model is well char- \ t \ 2 L t \ 2 

acterized by <J q =o. 5 (t) ~ f A , which is the same time dependency for | J \ \ J \ 
q = 0.5 of the PME [Eq. (8)]. 



Analytical derivation of anomalous diffusion in the associated 
Markov model of feed-forward DTQW. The relationship between 
our model and the PME can be explored using the decomposition 
method of Romanelli et al 40 ' 49 , in which the unitary evolution of a 
DTQW model is decomposed into Markovian and interference 



\b t+1 \ 

n +1 \ 



+ i 



-2a 1 



'+U 2 |^| 2 +2 A /i- 



(9) 
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Figure 2 | Anomalous slow diffusion of the feed- forward DTQW model. Its probability distribution at t = 10 7 step displayed in Panel (A) with running 
averaged over 10 data sets (light blue line) is fitted by the ^-Gaussian (7) with q = 0.5 (red line) to obtain the ^-generalized standard deviation o q {t) in 
Panel (B). Panel (C) shows the long-time evolution of the ^-generalized standard deviation a q (t) (green dots), which is well fitted by a q=Qm5 {t) ~ f A 
(red line). 
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where the two terms including jSj = ^Re ^gja* (b^j J are interference 

terms, and ^Re[z] is the real part of a complex number z. 

Neglecting the interference terms and introducing the abbrevia- 

|2 . | ,|2 

we get the associated Markovian 



tions V- = | 
model; 



and R*= m 



-L t+1 - 
L j-i ■ 



}(* 



(10) 



11 



The numerical simulation of the associated Markovian model is 
performed under initial conditions of (i^^o) = (^1^1) = l/4> anc ^ 
the typical probability distribution shown in Fig. 3 (A) is well fitted by 
the ^-Gaussian with q = 0.0. Furthermore, Fig. 3 (B) shows that the 
time evolution of o q {t) of the associated Markovian model is well 
fitted to o- q=0 0 (t) ~ f 33 , which again is the same time dependency as 
the PME for q = 0. 

It is known that the classical Markovian model, i.e. one without the 
interference terms of the homogeneous DTQW, satisfies the stand- 
ard heat equation in the continuous limit. Consequently, the assoc- 
iated asymptotic probability distribution is a standard Gaussian. This 
implies that the ballistic transport property of the homogeneous 
DTQW comes from the interference term 40 . We thus consider the 
continuous limit 50 of the associated Markovian model. 

We introduce the density p(x, t) and current j(x, t) as 



p(x,t)-- 



Kx,t) = (sj- 



f Ax, 



(12) 



where Ax is the difference of the nearest-neighbor sites. Taking a 
Taylor expansion of Eq. (10), we get 

a2 



8t 



p(x,t) 



1 8 , ^ 



= 0, 



in the diffusion limit, i.e., the quantity (Ax) 2 / At remains constant (set 
to unity here for simplicity) as At, Ax —> 0 with the one-step time 
difference At. In a similar manner, by expanding Eq. (11) and taking 
the diffusion limit, we obtain 



1 



2(\-p(x,t))dx 



P{x,t), 



(14) 



which implies a breakdown in Fick's first law (j oc —dpldx) and is 
the hallmark of anomalous diffusion. By substituting Eq. (14) into 
Eq. (13), we obtain the following nonlinear PDE: 



d , . 1 

— p[x,t) = ~ 

dt 2(1- p(x,t)) 2 



(15) 



Evaluating the asymptotic solution of this nonlinear PDE, after a 
long-time evolution, p(x, t) becomes much less than unity. As the 
rough approximation in this long-time limit, we have 1 — p ~ 1 and 
p 2 ~ 0, and Eq. (15) is thus well approximated by 



d / x 1 dl 2/ x 



(16) 



which is nothing but the PME in Eq. (6) with m = 2 (q = 0). We thus 
conclude that the approximated asymptotic solution of Eq. (15) is a 
^-Gaussian with q = 0. In addition, we can show that this result is 
mathematically valid by applying the asymptotic Lie symmetry 
method 51 (see Method). This method can give an equivalence 
between the asymptotic solution of the PDE and the analytically- 
solved one of the other PDE without analytically solving this PDE. 
Therefore, the associated Markovian model exhibits anomalous dif- 
fusion described by the PME in Eq. (6) with m = 2. This implies that 
the interference term of our model leads to the speed-up of the 



quantum walker o q = Q3 
Markovian model a, 



q = 0 ' 



f A compared to the associated 
and makes the zig-zag shape around 



the q- Gaussian distribution. 

In summary, we have proposed a feed-forward DTQW model Eq. 
(4) in which the coin operator depends on the coin states of the 
nearest-neighbor sites. We show that this model is experimentally 
feasible. Our feed-forward DTQW model asymptotically satisfies the 
PME for m = 1.5 (q = 0.5) and exhibits anomalous slow diffusion 



(13) a q=03 (t) ~ t 0A from the probability distribution and the time 



dependency of the standard deviation defined in the g- Gaussian 
distribution. 



Discussion 

In this section, we show that our results after the long-time numerical 
simulations have no initial coin dependence, and that the inter- 
ference term can be taken as the noise source in addition to the 
PME. First, while the above analysis uses the only fixed initial coin 
states as (flo»^o) = { a °i^°i) = (1/^^/2), we numerically confirm that 
there is almost no dependence of the initial coin state except for the 
trivial cases as follows. We have performed the several numerical 
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-400 
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Figure 3 | Anomalous slow diffusion of the associated Markovian model for the nonlinear quantum walk. Panel (A) shows the probability distribution 
of the associated Markovian model at t= 10 7 step (green dots) fitted by the ^-Gaussian, yielding q = 0.0and<Tg =o .o(f) = 283 (red line). Panel (B) shows the 
long-time evolution of the ^-generalized standard deviation a q (t) of the associated Markovian model (blue dots). It is well fitted by G q =o,o(t) ~ f 33 
(red line). 
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simulations for the initial state specified by (flgj^o 

(cos fin j v^sin fin j \/2) and = ^cosyjt j \fl,smyn j \fl 

with the real-parameter p and y ranging from 0 to 1. Note that the 
trivial cases, /? = 0.5, y = 0 and ft = 0, y = 0.5, lead to the localization 
of the probability distribution for any time, and we cannot define the 
parameter q for the trivial initial states. Figure 4 shows the numerical 
evaluation of the parameter q of ^-Gaussian distribution from the 
data at the two different time steps t = 10 6 and t = 10 7 , under the 
assumption to satisfy the stationary solution of the PME [Eqs. (7) and 
(8)]. The evaluated g -parameters for the various initial states are 
q = O.5+0 J47 except for the trivial cases. Therefore, we can conclude 
that our nonlinear model shows the anomalous slow diffusion to 
satisfy the PME with m~1.5 (q~0.5) without the initial state 
dependence. 

Finally, let us consider the difference between the probability dis- 
tribution of our model and the ^-Gaussian distribution with q = 0.5, 
as shown in Fig. 2 (B); the power spectrum of this difference exhibits 
a white noise as shown in Fig. 5. This power spectrum divided by the 
physical time scale t 0A may remain finite in the asymptotic case, 
which suggests that our nonlinear model may be mapped to the 
stochastic PME, i.e. the PME plus a white noise term, in the continu- 
ous limit. This stochasticity must come from the interference term. 
The problem of extracting the stochasticity from a deterministic 
process has been discussed in another context, that of Mori's noise 52 . 
Further analysis of this model may reveal the origin of the stochas- 
ticity. This is interesting as a purely mathematical problem of a 
stochastic nonlinear partial differential equation and for showing 
the relationship between the discrete model and its continuous limit. 

Methods 

In what follows, the solution of Eq. (15) is asymptotically identical to the solution of 
Eq. (16). This is mathematically equivalent to showing that the probability distri- 
bution 



„(*=<>) 



(*) 



Z((T q = o) 



u q = 0 



(17) 



is invariant under an asymptotic Lie symmetry 51 of the nonlinear partial differential 
equation (15). In other words, 



dtP = 



2(1- pf [2 



dxxP ~P dxxP 



(18) 



In Eq. (17), Z(a q=0 ) = 4a q=0 /3 is the normalization factor, and in what follows, the 
argument of this function is omitted where possible and d t p is denoted as p t for 
simplicity. 

We follow the asymptotic Lie symmetry method and notations in Ref. 5 1 . Under an 
infinitesimal transformation with the generator 




Figure 4 | Initial coin state dependence. Changing the parameters /3 and y, 
we numerically evaluate the parameter q of ^-Gaussian distribution for 225 

(cos fin j y/l, sin [in / y/l) 

and (a?,*; 



different initial states expressed by \ a 0 ,i 

cos yn j \/2,sin y% j y/lj . Note that the trivial cases, ft = 
0.5, y = 0 and /? = 0.5, y = 1, are not plotted. Our fitting result except for 
the trivial cases is q = 0.5^ 



0.0003 




4000 
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Figure 5 | The difference between the nonlinear model and the fit. The 

power spectrum of the difference between the probability distribution of 
our model and the ^-Gaussian with q = 0.5 at 10 7 step. To remove the 
effects of the expectation value, we replace x with x — 36.91 in the q- 
Gaussian with q = 0.5 [Eq. (7)]. 



X = Zd x + td t + <pd pi (19) 

that is 

Xh^-X = X + £^(x,t,p), 

t^t = t + £T(x,t,p), (20) 
p\-^f) = p + £(p(x,t,p), 

the function p(x, t) is mapped to a new function p(x,t), with 

p(x,t) = p(x,t)+e[(p- pj- p t T] p=p{xt) . (21) 

By applying this to the probability distribution Eq. (17), we see that the trans- 
formation X with £ = — x leaves Eq. (17) invariant if and only if 



(p = pj= -p x x = 



2x 2 



(22) 



Note that x = rj • t remains unrestricted at this stage because p (q 0) (x) does not 
explicitly depend on time t. Conversely, the function p(x) is invariant under 
X = — xd x + xd t + 2x 2 j (Zo 2 q =Q S jd p for any x if and only if p(x) is of the form given in 
Eq. (17). 

Following the general procedure for a Lie group analysis of differential equations 53 , 
the second prolongation of X is described by 

Y = X-\- x ¥ t d pt -\- x i , x dp X -\- x i , xx dp Xx . (23) 

The coefficients ^F*, and are defined as follows. Under an infinitesimal 
transformation of X, the partial derivatives are transformed as p x \-^p x -\-e x ¥ x , 
Pt^Pt + eWf, and p xx ^p xx + e x ¥ xx . We then readily obtain 



4x 

"z^To' 



ZoW 



<p P =o, 



=0. 



(24) 



The coefficients *F X , and y ¥ xx are then obtained by applying the prolongation 
formula (2.39) from Ref. 52: 



(<Pp-^pt=-npu 



x f x = < Pp + ( < Pp-£p)pp = 



4x 



Zo: 



+ Px, 



q = 0 



f 0.116 
-0.047- 



V* = (p xx + 2<p xpPx + (p pp p 2 x + (<p p -2^) p x: 

We note that Eq. (18) can be written as 

Pt = Ci{p x ) 2 + C 2 p xx 

with 



: 2(l-p) 2 ' Cl ~ 2(l-p)- 



Zo 



+ 2 Px3 



q = 0 



(25) 



(26) 



(27) 



(28) 



(29) 
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The asymptotic Lie symmetry condition 

Y(p t - d ( Px ) 2 - C 2Pxx ) = W f - 2C lPx W x - C 2 V XX - C[ <p(p x ) 2 - C' 2 (pp xx = 0 (30) 
with 



C, =d D Ci 



C' 2 =d 0 C 2 



1 



2(1-,)" yiZtf 



can be written in the following compact form: 

A 0 (x,t,p) +Ai (x,t,p)p x + A 2 (x,t,p) (p x f +A 2 (x,t,p)p„ = 0. 



(31) 



(32) 



When the condition in Eq. (30) is fulfilled, each A k (k = 0, 1, 2, 3) function must vanish 
separately in the asymptotic limit 



|p(*,f)|->0 for |x|->oo, 



(33) 



implying that the variance a q=0 also becomes infinity in the asymptotic limit from Eq. 
(17); 



G q = o->oo for |jc| —►OO- 



The function A 3 can be expressed as 

which must be nonzero as a q=0 — > oo, unless we choose 

rj= -2. 

Making this choice, X becomes 

2x 2 

X=td x -2td t + T -^d 0 



* ) 

l-p)Za 2 q _ 0 j' 



and A 3 reduces to 



(i-p)^ =0 

Thus, A 3 — > 0 as a q=0 — > oo. 

In a similar manner, A 0 , A ly and A 2 are given by 

2 p t _ 4 



{\-p)Zal_ 0 ' 



{1-pfZa* 



(l-pYZal =0 



(34) 



(35) 



(36) 



(37) 



(38) 



(39) 



and all become zero as (j q=0 — > oo. Therefore, we conclude that the distribution in Eq. 
(17) is an invariant solution for the transformation X of Eq. (37), which is an 
asymptotic symmetry for large \x\ of the nonlinear partial differential equation Eq. (18). 
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